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RESUME 

Le  champ  de  vitesse  du  sillage  d’une  sphere  hypersonique  est 
determine  a  I 'aide  d'un  double  peigne  transverse  de  sondes  ioniques 
pouvant  comporter  jusqu'a  huit  paires  de  sonde^.  Les  mesures  ont  porte 
sur  des  sillages  d'une  sphere  en  titane  de  2.7  pouces  de  diametre,  lancee 
a  14  500  pieds  a  la  seconde,  dans  un  corridor  de  tir  rempli  d'azot  . 

Deux  pressions  furent  utilisees  soit  7.6  et  20  torr.  La  vitesse  de 
convection  mesuree  se  compare  harmonieusement  avec  la  vitesse  moyenne  de 
sillage  deduite  de  1 'experience  des  etincelles  conslcutives  pour  de  I'air 
dans  des  conditions  similaires.  De  plus,  la  technique  des  peignes  de 
sondes  a  permis  d'etendre  la  determination  du  champ  de  vitesse  au  voi- 
sinage  immediat  du  projectile,  ce  qu'il  n'est  pas  possible  de  faire  a 
I'aide  de  la  technique  des  etincelles  consecutives  par  suite  de  la 
difficulte  de  former  des  etincelles  distinctes  dans  cette  region  du 
proche  sillage  en  raison  du  fort  degre  d 'ionisation  qui  y  regne. 

En  comparant  les  donnees  obtenues  a  7.6  et  a  20  torr,  on 
s'aper^oit  que  la  valeur  de  la  vitesse  de  convection  le  long  de  I'axe  de 
vol  augmente  avec  la  pression.  D'un  autre  cote,  si  I'on  considere  la 
largeur  de  sillage  en  vitesse  on  constate  qu'elle  est  plus  grande  a 
basse  pression,  du  moins  dans  le  proche  sillage.  Au  dela  de  300  ou  400 
diametres  de  projectile,  cependant,  les  donnees  aux  deux  pressions 
etudiees  ont  tendance  a  se  superposer.  Ces  observations  semblent  cor- 
roborer  les  informations  provenant  d'autres  sources  a  I'effet  que  le 
sillap,  a  7.6  torr  de  pression,  d'une  sphere  de  2.7  pouces,  n'est  pas 
complgtement  turbulent  avant  300  ou  400  diametres  de  projectile. 
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ABSTRACT 


The  technique  of  measuring  convection  velocity  in  the 
hypersonic  wake  with  a  pair  of  in-line  probes  has  been  applied  on  a  large 
scale  to  the  mapping  of  the  velocity  field  in  the  hypersonic  wake  through 
the  use  of  a  transverse  survey  array  containing  up  to  8  ion-probe  pairs. 
Measurements  are  reported  of  the  velocity  field  in  the  wakes  of  2.7  inch 
diameter  spheres  flown  at  14,500  feet/second  in  ballistic  range  atmos¬ 
pheres  at  7.6  torr  and  at  20  torr  of  nitrogen.  The  array  technique  leads 
to  convection  velocity  results  which  are  in  excellent  agreement  with  the 
mean  wake  velocity  data  obtained  by  the  sequential  spark  technique  in  an 
air  atmosphere.  In  addition,  the  technique  has  permitted  an  extension  of 
the  mapping  of  the  velocity  field  of  15,000  feet/second  spheres  to 
considerably  smaller  axial  distances  than  was  possible  with  the  spark 
method  because  of  the  difficulty- to  forming  distinct  sparks  at  the  higher 
levels  of  ambient  ionization  encountered  in  the  near  wake. 

A  comparison  of  the  data  obtained  at  7.6  torr  and  at  20 
torr  shows  that  the  amplitude  of  the  velocity  distribution  in  the  wake, 
as  defined  by  the  velocity  on  the  wake  axis,  is  higher  at  higher  pressure. 
Considering  the  wake  width  however,  it  is  found  that  the  width  of  the 
velocity  distribution  is  larger  at  the  lower  pressure  of  7.6  torr  than  is- 
the  case  at  20  torr,  at  least  in  the  near  wake.  However,  at  axial  dis¬ 
tances  greater  than  300  or  400  diameters,  the  data  for  the  two  pressures 
tend  to  overlap.  These  observations  appear  to  be  related  to  other 
information  indicating  that  for  2.7  inch  diameter  spheres  launched  at 
7.6  torr,  the  wake  is  not  fully  turbulent  until  about  300  or  400  diameters. 
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1.0  INTRODUCTION 


The  measurement  in  the  turbulent  wake  of  hypersonic  pro¬ 
jectiles  of  the  axial  and  radial  variation  of  the  velocity  distribution 
is  of  considerable  interest  for  the  development  of  turbulent  mixing 
models  for  the  wake.  Such  models  of  mixing,  in  combination  with  tempera¬ 
ture  data  and  reacting  chemistry  schemes,  are  required  to  predict  the 
persistence  of  various  chemical  species  in  the  wake  (References  1,  2). 

Recently  systematic  production  of  statistically  significant 
quantities  of  spatially-resolved  mean  wake  velocity  data  has  been  achieved 
using  the  sequential  Lpark  technique  (References  3,  4).  In  the  spark 
technique,  a  column  of  ionized  gas  is  created  along  an  arbitrary  path 
across  a  wake,  and  the  distortion  of  this  path  is  subsequently  recorded 
at  known  time  intervals  so  that  a  displacement  history  and  consequently 
a  mean  velocity  estimate  can  be  obtained.  Both  sphere  and  cone  wakes 
have  been  studied  over  a  wide  range  of  ambient  pressure  and  projectile 
velocity  including  both  laminar  and  turbulent  wake  flows.  One  of  the 
few  limitations  on  the  technique  occurs  when  the  ambient  density  is  so  low 
or  the  ionization  in  the  wake  is  so  high  that  an  electrical  discharge 
through  the  wake  produces  only  a  glow,  rather  than  a  well-defined  spark 
path.  For  example,  this  limitation  makes  it  impossible  to  measure  a 
velocity  profile  across  a  sphere  wake  at  axial  distances  less  than  about 
300  body  diameters  for  a  15,000  feet/second  hypersonic  sphere. 

Another  velocity-measuring  technique  can  be  devised  for  the 
ionized  hypersonic  wake  by  locating  two  electrostatic  probes  in  the  wake 
so  that  one  is  exactly  downstream  (or  axially  displaced)  from  the  other, 
and  measuring  the  time  it  takes  for  an  ionized  blob  (or  eddy)  of  fluid  to 
travel  or  be  convected  the  known  distance  between  the  two  probes.  As 
concerns  their  general  applicability,  such  techniques  employing  two-or- 
more  probes  in  axial  configuration  depend  on  the  wake  being  turbulent; 
while  they  cannot  produce  velocity  data  at  radial  distances  less  than  a 
body  radius,  they  do  permit  velocity  measurements  just  behind  the  pro¬ 
jectile.  Of  course,  being  based  on  correlation  methods,  such  techniques 
generally  demand  considerable  computer  effort.  With  a  single  pair  of 
probes,  data  can  only  be  obtained  at  one  arbitrary  value  of  radial  distance 
on  a  given  round,  and  consequently  a  large  number  of  rounds  are  required 
to  obtain  reasonable  statistics  at  any  particular  value  of  radial  distance. 
Up  to  the  present,  rather  liii:ited  quantities  of  convection  velocity  data 
have  been  reported  from  measurements  with  axial  arrays  of  up  to  5 
electrostatic  probes  on  large  hypersonic  sphere  wakes  (References  5,  6) 
and  with  pairs  of  hot  wires  or  cooled  film  anemometers  (References  7,  8) 
on  supersonic  sphere  wakes . 

The  sequential  spark  technique  is  believed  to  measure  a 
mean  wake  velocity.  There  is  some  reason  to  believe  that  wake  convection 
velocity  estimates  may  be  different  from  the  so-called  mean  wake  velocity 
measurements  obtained  at  the  same  location  in  the  wake  (References  8,  9, 

10) .  However,  the  wake  convection  velocity  data  measured  with  pairs  of* 
probes  have  so  far  been  found  to  be  in  general  agreement  with  other 
wake  velocity  data. 
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efficiency  of  the  two-probe  wake  velocity  measuring 

distributed 

across  the  wake  in  a  transverse  survey  array  configuration.  In  the 
present  paper  we  report  on  velocity  measurements  made  in  the  wakes  of 
diameter  spheres  flown  at  14.500  feet/second  in  the  ballistic 
range  facilities  at  the  Defence  Research  Establishment  Valcartier  (DREV) 
u.ing  transverse  survey  arrays  of  electrostatic  ion  probe  pairs. 


2.1  EXPERIMENTAL  TECHNIQUE 

^  iofi  probe  pair  consists  of  two  separate  cylindrical - 

hSghrfF^u^ri/"  ®  support  about  12  inches  in 

S^t^ofTL  mill  cylindrical  probe  consists  of  a  2  millimeter  long 

iS  and  f  !"■  g°ld  wire,  formed  by  stripping  the  outer 

jacket  and  -cflon  insulator  from  miniature  solid- jacketed  coaxial  cable  to 
expose  the  center  conductor.  Probes  are  normally  biased  to  collLt  iJns 
a  compensating  current  is  collected  by  the  remaining  jacket  and  b^thr 

that°a^  "®®*  ®®cl*  2-probe  pair  is  mounted  in  the  range  so 

between  the  tips  of  the  two  probes  woufd  be 
parallel  to  the  axis  of  flight  of  the  projectile.  Up  to  eight  such 
elemental  2-probe  axial  arrays  are  mounted  together  to  form  one  so-called 
transverse  survey  array'  (Figure  2).  The  separation  between  the  collating 
""/"dividual  2-probe  element  is  0.95  centimeter;  as  Regards  tS^  ® 
tirfiist  nlbe®  5^®  ®®c°"d  probe  is  directly  in  the  fhadow  of 

of  wakfielocitv  ^  tightly,  erroneous  estimates 

ot  wake  velocity  are  obtained.  Systematic  investigation  of  the  effect  of 

probe  diameter  and  probe  spacing  on  the  wake  velocity  estate  obtained  with 
pairs  of  probes  have  been  carried  out  by  Ghosh  and  Richard  (Reference  111 
using  a  plasma  jet  facility.  For  probes  of  the  size  used  in  the  present  * 
experiment,  their  results  shows  that  underestimation  occurs  only  for  probe 
separations  of  less  than  0.7  centimeter.  Experiments  at  DRE?  usin^proS 

confirmed  the  existence  of  sever^^flS 
underestimation  effects  for  small  probe  separation.  No  such  effects  are 
believed  to  exist  at  the  normal  probe  separation  distance  of  0.95  centimeter 
because  of  two  observations .  The  first  of  these  is  based  on  thfresCuT 
Richard  which  indicate  that  even  for  probes  of  diameter  ^p  to 
®°""iderably  larger  than  those  employed  here,  there  is 
effect  as  long  as  the  separation  distance  is  at  least  one  c^ntime^e^ 

The  second  is  the  fact  of  the  exceptionally  close  agreement  of  the  present 
results  with  those  obtained  by  the  sequential  spark  experiment.  ^ 

1  *  V,  2-probe  elements  in  the  survey  array  are  separated 

laterally  by  multiples  of  1.25  inches.  Both  symmetrical  and  asXetri^l 
arrangements  of  probes  (with  respect  to  the  flight  axis)  have  S  SJoJed 
The  prob®  signals  are  fed  from  the  probe  tips  into  one  of  the  16  pre-^  ^ 
amplifier  box  beneath  or  above  the  survey  array.  The  shape  of  the  pre- 
amplifier  box  was  determined  by  considerations  of  minimization  of  tL 
strength  of  the  reflected  bow  shock  wave  system  from  structure  in  the 
hypersonic  range.  The  range  itself  was  the  10  foot  inside  diameter 
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hypersonic  Range  5  i'acility,  which  consists  of  a  425  foot  long  evacuable 
free-flight  facility  along  which  various  models  can  be  flown  (Reference  12) . 
The  actual  free  diameter  of  the  range  is  reduced  at  the  experimental 
station  to  slightly  less  than  5  feet  by  the  presence  of  a  smoothly-profiled 
fiberglas  wedge  treatment  installed  on  the  walls.  This  treatment  has  been 
found  experimentally  (Reference  13)  to  significantly  attenuate  the  outgoing 
bow  shock  system  originating  with  the  projectile.  If  the  range  walls  were 
not  treated,  the  shock  system  would  reflect  at  the  walls  and  then  collapse 
toward  the  axis  of  the  tank;  in  the  process,  significant  disturbance  of 
wake  measurements  would  be  caused  (Reference  14) .  The  launcher  itself  is 
a  light  gas  gun  which  routinely  launches  saboted  models  up  to  4  inches  in 
diameter  at  velocities  of  14,500  feet/second.  In  the  flight  chamber, 
model  attitude  and  position  are  routinely  determined  by  a  system  of  flash 
X-rays.  This  position  data  permits  the  determination  of  the  flight  axis 
of  the  projectile  with  respect  to  the  probes.  The  probes  are  optically 
surveyed  in  place  in  the  range  before  the  launching. 

The  electrostatic  survey  arrays  have  normally  been  employed 
using  saturation  ion  current  collection,  the  probes  being  biased  at  minus 
4  volts  with  respect  to  ground.  (The  plasma  potential  in  the  wake  will 
be  slightly  positive  with  respect  to  ground).  The  bias  is  supplied  by  an 
operational  type  current-to-voltage  preamplifier,  which  tends  to  neutralize 
probe  capacitive  effects  and  which  drives  the  eight  feet  or  so  of  coaxial 
cable  needed  to  transmit  the  probe  signals  to  the  exterior  of  the  range 
tank.  The  preamplifier  employs  a  feedback  resistor  of  200  kilohms  and 
has  a  response  which  is  flat  between  dc  and  200  kilohertz,  falling 
smoothly  off  to  -3  db  between  200  and  400  kilohertz  (Figure  3) .  The  bias 
voltage  is  subtracted  in  a  type  'O'  unit,  and  the  residual  probe  signal 
is  amplified  by  a  Hewlett  Packard  467A  dc  power  amplifier  and  delivered 
by  means  of  170  feet  of  terminated  92  ohm  cable  to  the  recording  room. 
Recording  is  carried  out  by  means  of  Tektronix  551  and  555  double  beam 
oscilloscopes  and  Wollensak  35  millimeter  Fastax  cameras,  which  provide 
an  effective  sweep  rate  of  about  200  microseconds /centimeter.  For  velocity 
measurements,  each  probe  of  a  pair  is  recorded  using  one  beam  of  the  same 
oscilloscope;  the  two  probe  signals  are  recorded  in  synchronization. 

Timing  marks  are  introduced  by  means  of  Z-modulation  every  millisecond, 
these  assist  in  synchronization  and  film  velocity  measurement.  Film 
velocity  variation  amounts  to  a  fraction  of  1%  between  adjacent  1  milli¬ 
second  signal  segments.  Figures  4a  and  4b  show  typical  recordings  of 
pairs  of  signals  obtained  with  a  survey  array  at  pressures  of  7.6  and 
20  torr. 


The  convection-velocity-measuring  technique  using  a  pair 
of  electrostatic  ion  probes  to  detect  convected  ionization  fluctuations 
in  the  turbulent  hypersonic  wake  has  a  rather  large  advantage  compared 
to  other  applications  of  continuum- type  electrostatic  probes.  As  long  as 
the  probes  detect  a  signal  which  is  roughly  proportional  to  the  local 
ionization  being  swept  past  the  probes,  it  is  possible  to  neglect  the 
problem  of  exactly  interpreting  the  probe  current  in  terms  of  charge 
density  (Reference  15) .  In  effect  one  is  only  interested  in  the  time 
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delay  occuring  from  the  moment  a  given  signal  feature  is  seen  on  the 
upstream  probe  to  the  moment  the  same  feature  appears  on  the  downstream 
probe.  Of  course,  there  may  be  many  such  features  in  a  short  signal 
segment,  and  one  does  not  attempt  to  measure  an  average  time  delay 
directly.  Instead  the  two  signals  from  the  probes  are  cross-correlated, 
and  the  resulting  cross-correlation  curve  is  used  to  provide  the  desired 
estimate  of  wake  velocity  (Reference  12) . 


3.0  ANALYSIS 

The  first  stage  in  the  analysis  process  is  the  digitalization 
of  the  analog  signal  traces  recorded  in  pairs  on  35  millimeter  film.  To 
read  the  film  a  'positive*  must  be  made,  where  the  signal  appears  as  a 
transparent  trace  on  a  black  background.  This  positive  is  fed  into  an 
automatic  or  computer-controlled  film  reader  assembled  at  DREV  (Reference  16). 
In  the  reader  a  small  spot  of  light  is  displayed  on  a  precision  CRT  and 
then  imaged  onto  the  film.  Collection  optics  behind  the  film  directs  the 
transmitted  light  to  a  photomultiplier  tube.  At  the  same  time  some  light 
from  the  spot  on  the  CRT  is  removed  with  a  beam  splitter,  passed  through 
similar  optics  and  uniform  density  filter  and  measured  with  a  second  photo¬ 
multiplier.  The  output  signals  of  the  two  photomultipliers  are  compared, 
and  the  results  of  the  comparison  indicate  whether  the  density  of  the 
particular  spot  on  the  film  is  more  or  less  than  the  reference  density. 

The  results  of  the  comparison  are  fed  to  a  digital  computer  which  controls 
the  coordinates  of  the  spot  displayed  and  which  is  programmed  so  as  to 
conditionally  scan  the  film  to  locate  and  follow  the  signal  trace.  The 
digitalization  of  the  signal  trace  is  simply  a  consequence  of  the  signal 
tracing  process  of  the  reader.  The  data  is  outputed  on  magnetic  tape  and 
subsequently  printed  out  on  a  line  printer  by  the  XDS  computer,  at  which 
stage  any  anomalies  are  spotted  and  control  cards  for  analysis  instructions 
are  prepared. 

The  timing  marks  on  the  signal  traces  form  a  convenient 
common  time  reference  for  all  the  signal  pairs  recorded  on  a  given  round. 

The  duration  of  the  timing  marks  (about  10  microseconds)  is  short  compared 
to  the  frequency  content  of  the  probe  signals;  the  reader  simply  inter¬ 
polates  across  the  gap  in  the  trace  but  flags  the  points  at  the  same  time. 

The  reading  frequency  normally  employed  is  about  450-500  kilohertz,  so  that 
the  4  or  5  consecutive  flagged  points  appearing  every  500  points  are 
readily  recognized  as  timing  marks.  By  means  of  trigger,  delay  units,  and 
calibrated  high-sweep-speed  oscilloscopic  recording,  it  is  possible  to 
relate  certain  features  on  the  film  traces  to  the  time  at  which  the  pro¬ 
jectile  passed  the  measuring  station  and  consequently  fix  the  tim?  of 
occurrence  of  the  reference  timing  marks  with  respect  to  the  time  of 
passage  of  the  projectile.  The  timing  marks  also  permit  the  synchroni¬ 
zation  in  time  of  the  signals  recorded  from  a  pair  of  probes. 
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Several  considerations  enter  into  the  choice  of  the  length 
(in  time)  of  the  signal  samples  which  are  cross-correlated  during  the 
velocity  estimating  process  (Figures  4,  5).  Decay  of  the  charge  density 
in  the  wake  and  slowdown  in  the  wake  velocity  as  functions  of  time  or, 
equivalently,  distance  behind  the  projectile  introduce  respectively 
nonstationarity  in  the  amplitude  and  in  the  frequency  content  of  the 
random  current  signals  detected  by  the  electrostatic  probes.  The  signals 
can  be  rendered  quasi-stationary  by  considering  sufficiently  short  segments 
of  signal.  On  the  one  hand,  it  is  preferable  to  have  a  short  signal  in 
order  to  more  precisely  localize  the  region  of  the  wake  where  the  velocity 
estimate  is  made;  on  the  other  hand  the  result  is  more  precisely  estimated 
from  a  statistical  point  of  view  if  the  length  of  the  samples  which  are 
being  cross-correlated  can  be  increased.  With  large  piojectiles  such  as 
the  2.7  inch  diameter  spheres  being  fired  at  DREV,  the  variation  of  probe 
signal  amplitude  and  wake  velocity  are  considerably  slower  in  time  than 
that  experienced  with  smaller  models,  and  it  has  been  found  convenient  to 
use  signal  segment  lengths  of  0.5  millisecond  duration.  The  segments  are 
chosen  to  begin  or  to  end  on  the  timing  marks;  there  are  two  individual 
segments  defined  between  any  pair  of  consecutive  timing  marks.  For  2  7 
inch  diameter  spheres  flying  at  14,500  feet/second,  the  projectile  recedes 
from  the  measuring  station  at  about  66  diameters/millisecond,  equivalently 
the  effective  axial  distance  behind  the  projectile  over  which  an  individual 
0.5  millisecond  signal  segment  is  correlated  varies  by  about  33  body 
diameters.  At  axial  distances  greater  than  about  100  diameters  this  axial 
distance  variation  is  barely  significant.  The  velocity  estimate  is 
attributed  to  the  mid-point  of  the  signal  segment  or,  more  exactly,  to  the 
axial  distance  behind  the  projectile  corresponding  to  the  value  of  time  at 
this  point.  For  axial  distances  at  which  the  normalized  local  wake  velocity 
to  projectile  velocity  ratio  has  dropped  to  0.1,  the  actual  length  of  the 
(0.5  millisecond)  sample  of  wake  convected  past  the  probes  and  over  which 
the  cross-correlation  function  is  evaluated  is  about  3  diameters;  corre¬ 
spondingly,  when  the  ratio  drops  to  0.01,  the  length  of  wake  involved  is 
only  about  half  a  diameter. 

Theoretically  one  would  argue  that  the  convection  velocity 
measured  by  a  pair  of  probes  is  an  'averaged*  velocity.  In  practice, 
because  one  feature  in  a  probe  signal  may  dominate  other  features  and 
consequently  the  cross-correlation  function,  the  velocity  estimates 
obtained  by  such  a  pair  of  probes  frequently  tend  to  be  more  like  instan¬ 
taneous  estimates  of  wake  velocity  parallel  to  the  flight  axis,  rather 
than  averap  values.  It  is  well  to  point  out  here  that  measurements  made 
in  a  wake  by  a  stationary  probe  in  the  range  as  the  wake  flows  past  and 
over  the  probe  are  transient -type  measurements.  Tlie  whole  measurement  is 
opr  in  about  10  to  30  milliseconds.  This  is  in  considerable  contrast  to 
pe  case  of  the  study  of  a  steady-state  jet,  where  the  probes  can  be 
inserted  at  a  given  fixed  distance  from  the  throat  of  the  jet  and 
measurements  made  until  the  cross -correlation  function  is  as  well  averaged 
as  the  experimenter  requires.  Because  of  the  considerable  tendency  of  the 
prepnt  2-probe  velocity  estimates  towards  'instantaneous'  velocity 
estimates,  there  is  reapn  to  believe  that  the  rms  velocity  fluctuation 
data  derived  pom  the  distribution  of  velocity  data  points  obtained  by  the 
present  technique  may  be  considerably  closer  to  the  true  value  than  might 
otherwise  be  believed. 
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transverse  sS^ev  problems  were  encountered  with  the 

velocity  with  the  equi^alen^oH  ^PP^rent  when  measuring 

survey  array  of  ei^h^ri^^^  f  I"  effect,  the 

velocity  hi^ory  pfokler^o elements  allows  axial 

f  :i^,T  ir%  -V 


curve  P1P2,  gfneraLd^froftwr.^  cross- correlation  function 

correlated  ^Farh  Jfc  ^  ^  °  ^^Pled  probe  signals  which  have  been 

U  fall  “ 

interval  used  in  digitalization  of  i-h  integer  and  At  is  the  sampling 
time  lag  which  corresponds  to  tho  probe  signals.  Now  the  appropriate 
between  a  pai^  o^Ss  is  cha^aot  "  travel 

cross-corrLaiion  functio^P  S  K  "'">  '^^ich  the 

-  -- ™V5:  Jir 

SI  LtShi  u“f  "^^brcsiH• 

if  a,b  are  the  coordinates  of  a  comp^arpoinrof^S  °orra!a?lS  °thartS°"' 
line  passxng  fto™  the  point  (0,1)  tS  the  ?oint  (a, S,  Sarfsl^p^aSaJ  S 


M 


(b-l)a. 


If  |m|  is  plotted  as  a  function  of  T.  the  resultina  n„y„c,  .  -n  u 
minimum  at  the  point  of  tangencv  Rofono  J  curve  will  have  a 

the  program  was  thoroughly  toctoi  P“^  gs^ieral  application, 

manuLly  by  findlnr^S  poiJt  whfrr.R  ^“‘^"‘^/hich  had  already  been  analyzed 

the  computL  and  pfoS:d"c;Sa  taSSa^ei:  d°n  th^TS 
«re\‘han°'’5“‘"'?hSn‘"“  “  expariS:adSec\Mc-:n“aSV 

LXiralSe.oSt7:s?“S  SSi°Sa7d1s1nShMr™‘"^'j°"  “ 

curve  for  the  envelope  by  a  straight^ine  It  i 

if  the  error  introduLd  is  systemftic  oro  tMn  ^  5  evaluate 

certainty  is  that  this  eiro^irnegJ,- ^iM !  ^  ^ 

of  the  results  observed  because  o/thf  statistiS^^r"  scattering 

driven  fluctuations  in  the  velocity  LtimatL  """  turbulence- 

sHsriiilif SS  ==ir- 

been  analysed  separately.  n.e’'nSLS^'rL'r7nts  in”  till  LViS 
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considerable.  For  example,  the  fifteen  rounds  fired  at  7.6  torr  furnish, 
for  various  values  of  the  radial  distance  parameter  R/D  and  for  values 
of  the  axial  distance  X/D  behind  the  projectile  extending  from  0  to  about 
1500  diameters,  roughly  1200  independent  spatially-resolved  or  point 
esti'^ates  of  convection  velocity.  The  approach  adopted  in  order  to  deal 
with  this  data  consisted  of  dividing  the  results  into  sets  corresponding 
to  various  bands  of  axial  distances  behind  the  projectile.  The  data  in 
each  band  was  fitted  by  a  biparametric  gaussian  distribution  of  the  forr; 


V 

V 


V 
_o 

V 


R2 

r^ 


where  the  least  mean  squares  method  is  used  to  find  the  amplitude  and  the 
width  of  the  distribution.  These  two  parameters,  respectively  Vq/Voo  and 
r/D,  represent  the  ratios  of  the  axial  velocity  of  convection  to  the 
velocity  of  the  projectile,  and  of  the  convection  velocity  radius  of  the 
wake  to  the  diameter  of  the  projectile.  The  adequacy  of  the  gaussian 
distribution  for  such  curve  fitting  has  been  established  for  wake  velocity 
by  the  spark  technique  (References  4,  17).  One  of  the  practical  questions 
that  should  be  faced  in  fitting  this  distribution  is  that  of  deciding  what 
should  be  the  widths  of  the  bands  of  axial  distance  that  define  the  number 
of  data  points  in  a  given  set. 


Several  considerations  enter  into  the  choice  of  a  suitable 
width  for  the  axial  distance  bands  which  contain  the  sets  of  points  used 
to  define  the  gaussian  fits.  One  can  define  a  lower  limit  on  the  basis 
of  the  following  argument.  Each  point  measurement  of  velocity  obtained 
with  a  survey  array  has  associated  two  coordinates:  a  radial  distance 
R/D,  which  is  essentially  exact,  and  an  axial  distance  behind  the 
projectile  X/D,  which  is  imprecise  by  about  ±  16,  based  on  the  fact  that 
the  basic  correlation  interval  is  0.5  milliseconds  or  32  body  diameters. 

It  would  be  futile  to  attempt  to  use  bands  narrower  than  32  diameters  for 
the  simple  reason  that  the  narrower  the  band,  the  fewer  the  number  of  data 
points  in  the  set. 


It  is  important  to  realize  that  one  of  the  parameters  that 
is  established  by  the  fit  is  the  axial  value  of  the  convection  velocity. 

For  physical  reasons,  it  is  impossible  to  get  closer  to  the  wake  axis  with 
a  pairs  of  probes  than  a  radial  distance  of  0.5  diameter;  most  of  the  data 
corresponds  to  radial  distance  values  greater  than  0.8  diameter.  In  order 
to  increase  the  statistical  precision  of  the  gaussian  fit  to  the  velocity 
data,  it  is  desirable  to  utilize  the  largest  possible  set  of  data  points, 
without  sacrificing  too  much  the  axial  resolution  of  the  measurements. 

One  can  ask  the  question  as  to  whether  there  exists  some 
optimum  bandwidth  which  would  reconcile  both,  precision  and  resolution. 

The  question  has  been  tackled  experimentally  on  the  XDS  computer  using 
the  actual  experimental  data.  Figure  6a  shows  a  comparison  of  calculations 
of  axial  velocity  using  three  diff3rent  bandwidths .  In  the  first  case 
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diameters  with  no  overlap 

the  band  wL  Vi  second  case  (open  circles)  the  width  of 

a  to?al  wiSh  of  "  an  overlap  of  15  diameters  at  each  end  for 

a  total  width  of  60  diameters.  Finall)  in  the  third  and  last  case  fsolid 

As^'is^eiident^f^^^  diameters  with  no  overlip. 

nLitic^nr^w  different  bandwidths  give  ^ 

?he  iidth  nf  thus  indicating  that  within  certain  limits 

ieLc  as  a  WtJon  T  -  the  profile  of  convecti'^^ 

wit^riLS  to  Tfi  distance.  Similar  findings  were  made 

itn  regard  to  the  radius  of  the  velocity  wake  (Figure  6b) . 


4.0  RESULTS 

f-rnm  typical  radial  profile  of  convection  velocity  derived 

slhlrt  ‘^;^"'®ter  band  of  data  centered  at  420  diameters  ^behind  a 

p  re  IS  shown  in  Figure  7.  The  curve  passing  throuvh  ^he  data  is 
gaussi™  having  ^  amplltudn  on  the  wak^  axis\/v.  ?  O^S.  a^d  a  radius 

Circles  indicate  the  mean  values  obtained  by  averaging 
all  the  data  points  in  radial  sectors  of  width  0.2  diameter  take7cor 

thrialu^orth^n^^®  normalized  by 

i  f  ^  normalized  wake  velocity  on  the  axis  (Vq/Vco)  .  Figure  7 

certain  features  common  to  all  such  radial  velocity 

^  0"e  notes  that  the  velocity  estimates  are  not  uniformly  dis- 

ributed  as  a  function  of  R/D:  there  are  none  at  values  of  R/D  iLs  than 

0.8,  and  very  few  corresponding  to  radial  distances  greater  tLn  the 

=  2-28)  •  This  haf  tL  resuU  that  the 
wake  radius  is  not  very  well  determined,  particularly  at  larger  axial  dis 

h^'c  nl^fof  Ltf  JoinJ:\n^" 

the  viciniy  of  the  wake  axis  to  ensure  that  the  axial  value  of  the  assumed 
gaussian  distribution  of  velocity  is  determinable  with  acceptable  precision. 

V  /V  Quantitatively,  the  precision  of  the  derived  parameters 

Vo/V«  and  r/D  can  be  displayed  with  the  aid  of  contours  of  levels  of 

^  confidence  contours  encircling  a 

oaraLterfof  the^^  ^  coordinates  are  the  amplitude  and  width 

parameters  of  the  gaussian  curve  in  Figure  7.  The  'percentage'  attached 

valuef  0^"  probability  that  the  true 

values  of  Vq/Voo  and  r/D  are  enclosed  within  that  contour.  We  have  chosen 

somewhat  arbitrarily  to  use  confidence  contours  at  90%  in  order  to  specify 
the  precision  attached  to  the  results  of  the  radial  profile  curve-fitting 

'contour'  to  each  fit,  it  is  simpler  to  give  only 
V  projection  of  these  confidence  contours  on  the 

rectangle  °  dimensions  of  the  circumscribing 

rectangle.  These  limits  of  confidence  are  considered  to  be  relatively 

conservative.  For  example,  the  specification  of  a  confidence  level  of  90% 
is  equivalent,  in  the  case  of  a  gaussian  error  curve,  to  specifying  error 
limits  at  ±  1.64  a,  where  a  is  the  standard  deviation.  ® 
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a  2  7  behavior  of  the  velocity  on  the  axis  of  the  wake  behind 

nitroair\  feet/second  in  a  7.6  torr 

nitrogen  atmosphere,  previously  given  in  Figure  6,  is  plotted  aeain  in 

the  90%  confidence  limits  added  on  certain  points^ 

given;  these  correspond  to  the  previously- 
mentioned  imprecision  of  ±  16  diameters  in  the  value  of  the  axial  distance 
to  which  a  convection  velocity  estimate  belongs.  Figure  9a  shows  the 

watri°H  °^^he  velocity  profi!e  0?  the 

aS  ?Annendr''r?®°“'i^  behavior  of  the  velocity  on  the  wake 

Si^ani  S!  J  furnishes  tabulations  of  both  the  wake  velocity  on  the 
axis  and  the  wake  radius  as  a  function  of  axial  distance  from  2  7  inrh 
sphere  firings  at  7.6  torr,  for  several  bandwidths) . 

Fieure  9  present  data  similar  to  that  shown  in 

in  that  they  give  respectively  wake  velocity  on  the  axis  and  the 

ZlTU  ^  distance,  ex^epr^hat 

here  the  ambient  pressure  in  the  range  was  20  torr  nitrogen.  (The  data 
are  tabulated  in  Appendix  B) .  ^ 

.1-n.^c  w  ^^8ure  9  for  firings  in  7.6  torr  nitrogen 

r  for  firings  at  20  torr,  are  plotted  together 

alnna^th!  f^rst  the  behavior  of  the  wake  convection  velocity 

f  rfori  u°"®  profile  corresponding  to  ^ 

rP«J)  -  t  "  w  f°rr -inches)  falls  below  that  corresponding  to  20  torr 
L??  “  torr-inches).  The  difference  between  the  two  profiles  is 
systematic  and  statistically  significant,  and  the  direction  of  the  dif- 
erence  indicates  that  an  increase  in  ambient  pressure  entails  an  increase 

estimates  of  the  mean  wake  velocity  on  the  axis  determined  by  Lahaye  et 
al.  (Reference  4)  and  Lahaye  (Reference  17)  using  the  sequential  spark 
technique  (and  15,000  feet/second  spherical  projectiles  launched  in  air 

s  experimental  conditions 

similar  to  those  described  here  except  that  the  diameter  of  the  spherical 
projectiles  used  did  not  exceed  one  inch.  0D\iously,  the  spirk  rLults 
agree  very  well  with  those  obtained  using  arrays  of  probes,  despite  any 
difference  in  atmospheric  composition.  One  concludes  that  the  Lray 
measurements  confirm  the  spark  data,  and  that  there  is  no  difference 
between  the  probe  convection  velocity  and  the  mean  wake  velocity. 
Additionally,  the  probe  array  technique  has  permitted  a  mapping  of  the 
wake  velocity  field  In  the  region  of  axial  distance  less  thajfoo  dimeters 
behind  the  spheres.  As  is  indicated  by  the  straight  line  of  slope  ToJ 
he  graph,  the  decay  of  the  wake  convection  velocity  on  the  axis^follows 

di^eterX^*"  °  distances  X/D  between  300  and  1000 


Concerning  the  velocity  radius  of  the  wake,  as  has  been 
previously  discussed,  the  present  array  techniques  have  not  permitted  us 
to  determine  the  wake  radius  with  precision.  In  effect,  the^problem  is 
attributable  in  large  measure  to  a  distribution  of  the  measurLents  over 
radial  distance  which  does  not  lend  itself  to  wake  width  determination 
Relative  to  the  size  of  the  wake,  the  individual  velocity  estimates  Ze 
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at  large  axial  distances 

Let  thatLn  FiJ^fnhLh  recognized,  it  still  remains  the 

thP  ?  Figure  lib  there  appeals  to  be  a  difference  in  the  size  of 

the  wake  velocity  width  at  7.6  torr  and  the  size  of  th^wake  at  20  Lrr 

environment  in  which  a  hypersonic  wake  develops  is  charac- 
lo  the  description  of  flow  fields,  the  effects  of  variabirden^tro; 

of=h:^“^dSL;^r;L^r:L^^L^^rL^1;L^^ru^ 

pR/D  ^ 


(Rhd/»)^ 


pR/ 

Wo 


ydy. 


sysjOT,''“““'  laboratory  coordinate 

p  is  the  wake  density 

Poo  fa  the  density  at  free  stream 

y  is  a  dummy  variable  of  integration. 

tend  to  behave  in  the  san,e  fashion  as  IS  sp«d  llake  ?he  tositv 

carrle7oit1t'DREv"^  t°Se  iaLr^rfT^Jerd''  7““”^ 

14  i;nn  ^  .  waxes  ot  aj,7  inch  diameter  spheres  flown  at 

in’the  pro  e  dL°gs  o"f  a  N«o'spe”aUsf's 

flows  (Reference  19).  J’Pecialist  s  Conference  on  turbulent  shear 

oniv  Howarth-Dorodnitsyn  transformation  was  originally  annlied 

only  to  compressible  laminar  flows;  however.  Laufer  (Reference  20?^hao  1 
that  It  can  be  equally  applied  to  compressiMe  tuibuient  knur.  ?r 

^adiur^  rr  *“  thrHoi^th-DoStt  ^ 

nitvnLn  A  ^  diameter  sphere  at  14.500  feet/second  in  7  6  torr  ^ 

nitrogen.  As  one  can  see  the  modification  introduced  by  this  cooLiZL 
transforation  is  far  from  negligible.  For  example  at  an  axLrdlstance 
f  100  diameters,  a  velocity  estimate  obtained  at  R/D  =  2  is  displaced  to 
Rhd/D  -  1  after  the  transformation.  displaced  to 
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The  axial  profile  of  convection  velocity  on  the  wake  axis  in 
the  transformed  coordinate  system  as  well  as  the  transformed  wake  velocity 
width  are  shown  in  Figure  13.  These  curves  have  been  obtained  by  fitting 
the  transformed  radial  velocity  data  witli  a  least  mean  squares  routine 
exactly  as  was  previously  done  on  the  untransformed  data.  The  remarkable 
fact  concerning  the  two  profiles  of  Figure  13  is  that  they  apparently 
ollow  the  asymptotic  laws  for  turbulent  wakes,  i.e.  a  decay  of  -  2/3  for 
the  velocity  on  the  wake  axis  and  a  growth  law  of  +  1/3  for  the  wake 
width.  The  points  in  the  region  of  axial  distance  of  1000  diameters 
depart  from  the  general  trend,  howevei .  It  should  be  pointed  out  that 
these  latter  points  are  the  least  precise  because  of  the  small  number  of 
measurements  available  in  this  region. 

Unfortunately,  owing  to  the  unavailability  of  density  data 
or  wakes  at  20  torr  it  has  not  been  possible  to  apply  the  transformation 
to  the  convection  velocity  data  obtained  at  this  pressure. 


5.0  DISCUSSION 

A  comparison  of  the  experimental  measurements  presented  in  the 
preening  section  raises  several  questions  of  considerable  importance. 

It  will  be  recalled  that  the  results  concerning  wake  convection  velocities 
were  all  obtained  under  practically  the  same  experimental  conditions, 
except  for  pressure,  which  alternated  between  7.6  and  20  torr  nitrogen. 

In  all  cases,  the  projectile  consisted  of  a  hollow  titanium  sphere  of 
diameter  2.7  inches,  flown  at  a  velocity  of  approximately  14,500  feet/ 
second.  Figures  11a  and  11b  summarize  graphically  the  influence  of 
pressure  on  the  behavior  with  axial  distance  of  the  convection  velocity 
on  the  wake  axis  and  of  the  width  of  the  wake.  One  notes  in  the  latter 
case  a  somewhat  intriguing  behavior.  At  the  lower  pressure  (7.6  torr  N2I 
the  size  of  the  wake  is  relatively  constant  until  about  700  diameters, 
after  which  distance  the  wake  continues  to  grow  but  apparently  in  a 
isorganized  fashion.  Contrastingly,  at  higher  pressure  (20  torr  N2I, 
the  wake  width  begins  to  grow  in  the  neighborhood  of  200  diameters 
apparently  follow^ Mg  a  1/3  power  rate.  Another  curious  feature  is ’that 
or  axial  distances  less  than  300  or  400  diameters  behind  the  projectile, 
the  width  of  the  velocity  wake  at  20  torr  is  smaller  than  that  obtained 
at  /.6  torr.  One  can  ask  the  question  as  to  whether  these  differences 
are  due  to  some  qualitative  differences  between  the  turbulence  at  the  two 
pressures  or  if  possibly  the  differences  have  something  to  do  with  the 
laminar  turbulent  transition. 

th.t  schlieren  studies  (Reference  21)  it  appears 

hat  behind  blunt  bodies,  transition  occurs  abruptly  from  straight  laminar 
flow  to  well-developed  turbulence  within  a  few  axial  body  diameters.  In 
the  case  of  the  present  measurements  at  7.6  torr,  the  Reynolds  number 
based  on  sphere  diameter  is  about  2  x  10^  and  from  the  data  of  Wilson 
(Reference  21),  the  normalized  transition  data  should  be  about  7  diameters 
behind  the  projectile.  Thus  the  wake  should  be  turbulent  at  the  two 
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pressures  considered  in  this  paper.  Nevertheless  the  behavior  of  the 
signals  (Figure  4)  tends  to  raise  doubts  as  to  the  character  of  the 
^  exhibited  by  the  wake  at  7,6  torr,  at  least  in  the  first  200 

or  300  diameters.  In  the  near  wake  (X/D  <  100),  the  signals  recorded  at 
7.6  torr  contain  very  few  fluctuations  at  high  frequencies  compared  to 
e  torr  signals;  in  fact,  at  7.6  torr,  high  frequencies  make  their 
appearance  only  between  100  and  300  diameters.  This  observation  is 
corroborated  quantitatively  by  the  relatively  high  number  of  estimates 
o  very  large  space  scales  (larger  than  the  sphere  diameter)  obtained, 
in  the  near  wake  at  7.6  torr,  whereas  there  are  hardly  any  such  at 
20  torr  (Reference  22) . 


apparent,  then,  that  there  are  several  well-defined 
differences  in  the  observed  behavior  of  the  wake  at  7.6  torr  and  at  20 
torr.  One  could  tend  to  conclude  from  these  differences  that  the  wake 
at  7.6  torr  is  transitionary  in  the  region  near  to  the  projectile  and 
does  not  become  completely  turbulent  until  about  300  diameters,  thereby 
admitting  a  contradiction  between  the  present  results  and  those  of 
Wilson.  However,  examination  of  Figure  13  leads  us  to  reconsider  this 
tentative  conclusion.  In  effect,  if  one  suppresses  the  effect  of  com¬ 
pressibility  from  the  data  at  7.6  torr  by  applying  the  coordinate  trans¬ 
formation  of  Howarth-Dorodnitsyn,  one  obtains  an  axial  behavior  of  the 
convection  velocity  measured  on  the  wake  axis  and  of  the  wake  width  which 
seems  to  conform  to  the  asymptotic  laws  characteristic  of  low  speed 
turbulent  wakes.  Thus  on  the  basis  of  Figure  13,  the  wake  of  a  sphere 
tlown  in  a  nitrogen  atmosphere  at  7.6  torr  would  become  completely  tur¬ 
bulent  around  150  diameters  and  the  turbulent  flow  would  be  self-similar 
However  after  Corrsin  and  Kistler  (Reference  23) ,  a  minimum  condition 
for  self-similarity  is  that  the  ratio  of  the  standard  deviation  of  the 
instantaneous  position  of  the  turbulent  front  in  the  wake  to  the  mean 
position  must  be  constant.  Measurements  of  intermittency  at  DREV 
indicate  that  between  100  and  400  diameters,  this  ratio  can  be  represented 
Dy  a  straight  line  of  positive  slope  (Reference  24) . 

As  one  can  see,  the  interpretation  of  the  results  presented  in 
this  article  is  not  an  easy  task.  At  the  same  time  we  have  no  doubt  that 
these  results  conform  well  to  physical  reality. 


6.0  CONCLUSION 


The  technique  of  measuring  convection  velocity  in  the  hyper¬ 
sonic  wake  with  a  pair  of  in-line  probes  has  been  applied  on  a  large 
scale  to  the  mapping  of  the  velocity  field  in  thi:  hypersonic  wake  through 
the  use  of  a  transverse  survey  array  containing  up  to  8  ion-probe  pairs 
Measurements  are  reported  of  the  velocity  field  in  the  wakes  of  2.7  inch 
diameter  spheres  flown  at  14,500  feet/sccond  in  ballistic  range  atmospheres 
at  7.6  torr  and  at  20  torr  in  nitrogen.  The  array  technique  leads  to 
convection  velocity  results  which  are  in  excellent  agreement  with  the  mean 
wake  velocity  data  obtained  by  the  sequential  spark  technique  in  an  air 
atmosphere.  In  addition,  the  technique  has  permitted  an  extension  of  the 
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mapping  of  the  velocity  field  of  14,500  feet/second  spheres  to  con- 
smaller  axial  distances  than  was  possible  with  the  spark 
method  because  of  the  difficulty  of  forming  sparks  at  the  higher  levels 
ot  ambient  ionization  encountered  in  the  near  wake. 

A  comparison  of  the  data  obtained  at  7.6  torr  and  at  20 
torr  shows  that  the  amplitude  of  the  velocity  distribution  in  the  wake, 
as  defined  by  the  velocity  on  the  wake  axis,  is  higher  a*-  higher  pressure 
Considering  the  wake  width  however,  it  is  found  that  the  width  o/the 
velocity  distribution  is  larger  at  the  lower  pressure  of  7.6  torr  than 
IS  the  case  at  20  torr,  at  least  in  the  near  wake.  However,  at  axial 
distances  greater  than  300  or  400  diameters,  the  data  for  the  two  pressures 
tend  to  overlap.  These  and  other  observations  on  the  wakes  at  7.6^and  20 
torr  tend  to  indicate  that  the  wake  behind  a  2.7  inch  diameter  sphere  at 
7.6  torr  may  not  be  completely  turbulent  before  300  or  400  diameters 
behind  the  projectile. 
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APPENDIX  A 


This  appendix  contains  three  tables  giving  all  the  pertinent 
numerical  values  pertaining  to  the  fits  of  the  radial  profiles  of  the 
convection  velocity  at  7.6  torr  (20  torr-inches) .  Each  table  is  charac¬ 
terized  by  a  three  number  combination  (per  example  60-30-60)  describing 
the  way  the  axial  bands  are  built.  The  first  number  designates  the  axial 
coordinate  of  the  center  of  the  first  band,  the  second  one  the  distance 
between  the  center  of  two  consecutive  bands  and  the  third  one  the  width 
of  the  bands.  It  is  to  be  noted  that  the  bands  will  overlap  if  the  third 
number  is  greater  than  the  second  one.  The  three  tables  give  the  same 
quantities  in  the  same  order  that  is  from  left  to  right,  the  axial  dis¬ 
tance  (X/D),  the  fitted  value  of  the  convection  velocity  on  the  wake  axis 
(Vo/Voo) ,  the  fitted  value  of  the  velocity  wake  width  (r/D) ,  the  residual 
si^  of  squares  (R.S.S.),  the  standard  deviation  of  the  fit  (S.D.)  and 
finally  the  number  of  measurements  forming  the  radial  profile  (N) . 
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TABLE  A- I 

NUMERICAL  DATA  RELATED  TO  THE  FITS  AT  7.6  TORR 
(Combination:  60-30-30) 


Vq/V^  r/D  R.S.S. 


60.0 

0.144 

90.0 

0.103 

120.0 

0.110 

150.0 

0.090 

180.0 

0.095 

210.0 

0.082 

240.0 

0.  066 

270.0 

0.072 

300.0 

0.058 

330.0 

0.054 

360.0 

0.055 

390.0 

0.053 

420.0 

0.048 

450.0 

0.043 

480.0 

0.042 

510.0 

0.037 

540.0 

0.036 

570.0 

0.036 

600.0 

0.039 

630.0 

0.034 

660.0 

0.033 

690.0 

0.037 

720.0 

0.028 

750.0 

0.030 

780.0 

0.031 

810.  0 

0.025 

840.0 

0.022 

870.0 

0.024 

900,0 

■  0.019 

930.0 

0.020 

960.0 

0.018 

990.0 

0.023 

1020.0 

0.022 

1050.0 

0.019 

1080.0 

0.019 

1110.0 

0.017 

1140.0 

0.016 

2.01 

9. 8625F“3 

3.75 

2.2313F'‘2 

2.16 

8. 8263F“3 

2.68 

1. 3551P“2 

2.35 

1.0496F“2 

2.37 

8. 1595F“3 

2.83 

8.5724F“3 

2.10 

5. 5282F“3 

2.89 

1.1845F“2 

3.37 

1.0560F"2 

2.43 

5.9887F"3 

2.23 

6. 5648F"3 

2.35 

3.5921F“3 

2.42 

6.3150F"3 

2.67 

5. 5454F"3 

2.94 

4. 1791F~3 

2.  81 

3. 2356F~3 

2.51 

3.0384F”3 

2. 51 

8.2016F~3 

2.48 

1.5223F~3 

2.86 

1.8317F~3 

2. 33 

1.6360F~3 

3.54 

1. 48  76F~3 

2. 77 

4.6519F~4 

3.36 

5. 9006F~4 

3.84 

8.0075r”4 

3.21 

1.3622F~4 

3.84 

4. 9932F~4 

6. 55 

2. 9015F"4 

4.26 

2.a33ir“4 

4.13 

1.4467^"4 

3.05 

1. 8166F~4 

2.79 

1 . 0384F“4 

3.  94 

1. 3120F“4 

4. 64 

2. 3215r~5 

5.38 

1. 3693r“4 

5,62 

6.6355'’~5 

2.0272F“2 

26 

2. 6829F“2 

33 

1 . 5240F“2 

40 

1.8640F“2 

41 

1 . 6199F“2 

42 

1 . 2408F“2 

55 

1 . 2840F”2 

54 

1. 0026F”2 

57 

1 .4950F”2 

55 

1. 3611F"2 

59 

1 . 1055F"2 

51 

1..1026F"2 

56 

8 . 9345F“3 

47 

1 . 0814F“2 

56 

1 . 0638F~2 

51 

8. 7168F~3 

57 

8.6745F~3 

45 

8 . 5055F~3 

44 

1.381  1F“2 

45 

6. 0934F~3 

43 

7.339  aF“3 

6 

6. 5615F“3 

40 

7. 5640F“3 

28 

5. 2311F~3 

19 

6. 7372F~3 

15 

7. 5628F“3 

16 

3. 6908F“3 

12 

6. 7374F~3 

13 

4. 9172F“3 

14 

4. 4985F“3 

16 

3. 3359F“3 

15 

4.7652F“3 

10 

3. 3967F“3 

11 

3. 3065F-3 

14 

2. 7818F~3 

5 

4 . 4229F“3 

9 

4. 7030F“3 

5 
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TABLE  A-III 


NUMERICAL  DATA  RELATED  TO  THE  FITS  AT  7.6  TORR 


(Combinatjon:  100-100-100) 


X/D  Vq/Voo  r/D  R.S.S. 


100.0 

0.111 

200.0 

0.082 

300.0 

0.060 

400.0 

0.052 

500.0 

0.038 

600.0 

0.036 

700.0 

0.033 

800.0 

0.026 

900.0 

0.021 

1000.0 

0.021 

1100.0 

0.017 

2.40  6.6469£’“2 
2.47  3.7477r“2 
2.67  3.1353r“2 
2.28  1.9919F’“2 
2.87  1.4343F“2 
2.50  1.3040F“2 
2.74  4.9423F“3 
3.43  2.3607F”3 
4.24  1.2168F”3 
3.19  4.4170F“4 
5.02  3.0536F"4 


2.3342F“2  124 
1.5401F“2  160 
1.2983F“2  188 
1.0730F“2  175 
9.4978F“3  161 
9.9018F“3  135 
6,9953F“3  103 
6,5515F“3  57 
5,1999F“3  47 
3,5525F"3  37 
3.•)02^^F“3  30 
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APPENDIX  B 


This  appendix  contains  only  one  table  giving  all  the  pertinent 
numerical  values  pertaining  to  the  fits  of  the  radial  profiles  of  the 
convection  velocity  at  20  torr  (54  torr-inches) .  The  various  quantities 
and  symbols  used  in  this  appendix  are  described  in  appendix  A. 
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TABLE  B-I 


NUMERICAL  DATA  RELATED  TO  THE  FITS  AT  20  TORR 
(Combination:  60-30-60) 


X/D 

Vo/V„ 

r/D 

R.S.S. 

S.D. 

N 

60.0 

0.157 

1.87 

2. 4788'’"3 

1. 5744r"2 

12 

90.0 

0.168 

1.68 

4. 3626F“3 

1 . 476  SP"2 

22 

120.0 

0.163 

1 .63 

9. 0895r"3 

1 . 8348F”2 

29 

150.0 

0.147 

1.74 

1.3989r"2 

1.9992F"2 

37 

leo.o 

0.140 

1.73 

1.8414f’“2 

1 . 9002F”2 

53 

210.0 

0.115 

1.90 

2. 5483r'’2 

1.9358F"2 

70 

240.0 

0.106 

1 .84 

2. 2261F“2 

1.6894F'’2 

270.0 

0.111 

1 . 68 

1. 7392F“2 

1.4653F"2 

83 

300.0 

0.099 

1.75 

1. 7834F"2 

1 . 5219F“2 

79 

330.0 

0.087 

1.80 

1.9884F“2 

1.4947F"2 

91 

360.0 

0.078 

1.88 

1. 5048r"2 

1. 2931F“2 

92 

390.0 

0.070 

2.10 

9.0451F"3 

9.8620F"3 

95 

420.0 

0.062 

2.35 

1.3n03F"2 

1.1695F"2 

450.0 

0.054 

2.46 

1.3/13F“2 

1.1769F"2 

480.0 

0.049 

1.0039F“2 

1.0019r“2 

510.0 

0.049 

2.47 

8.7003F"3 

1.0177F“2 

86 

540.0 

0.052 

2.28 

8. 8248F”3 

1.0311F“2 

85 

570.0 

0.040 

2.65 

9. 5620F’'3 

1.0799F“2 

84 

600.0 

0.043 

2.91 

8. 2114F"3 

9.499  2F"3 

93 

630.0 

0.047 

2.36 

8. 3771F"3 

9.5423F“3 

94 

660.0 

0.046 

2.33 

1. 1595F"2 

1.1226F"2 

94 

690.0 

0.041 

2.68 

8.7578F"3 

1.1029F"2 

74 

720.0 

0.039 

2.91 

6. 3815F"3 

9. 7594F"3 

69 

750.0 

0.030 

2.95 

4. 9787F"3 

8. 7519F"3 

67 

780.0 

0.033 

3.34 

4.3985r"3 

8.1024r”3 

69 

810.0 

0.033 

3.08 

4. 2721F"3 

7.9262F"3 

840.0 

0.035 

2.82 

3. 3159F"3 

7. 6950F"3 

58 

870.0 

0.036 

2.64 

2.6702F"3 

7.4585F"3 

900.0 

0.035 

2.71 

1.8965F"3 

6. 2857F"3 

930.0 

0.032 

3.23 

1. 5572r"3 

5.6957F"3 

50 

960.0 

0.032 

3.50 

1. 2908P’"3 

6.0729F"3 

37 

990.0 

0.032 

3. 44 

1 . 4303F"3 

6.4859r“3 

36 

1020.0 

0.029 

3.62 

8.3188F“4 

5.4507F“3 

1050.0 

0.030 

3.45 

6. 9908F"4 

5. 5132F"3 

25 

1080.0 

0.033 

2.93 

5. 1690F"4 

4. 9613F“3 

23 

1110.0 

0.030 

3.27 

3. 7867F"4 

5. 2008F"3 

16 

1140.0 

0.027 

3.75 

3. 5112F‘4 

4.  f 846F"3 

18 

1170.0 

0.028 

3.86 

3. 9720F“4 

6.009ir"3 

13 

1200.0 

0.035 

2.6  4 

2. 2592r"4 

5.314  1^"3 
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